Dynamics of the decay r\ —> 37r 



S. Teige a , G. S. Adams 1 , T. Adams h '\ Z. Bar- Yarn 6 , 
J. M. Bishop h , V. A. Bodyagin f , D. S. Brown e N. M. Cason h , 

S. U. Chung c , J. P. Cummings e ' 2 , K. Danyo c , 
A. I. Demianov f , S. P. Denisov d , V. Dorofeev d , J. P. Dowd e , 
A. R. Dzierba a , P. Eugenio e ' 3 , J. Gunter a , R. Hackenburg c , 

E. I. Ivanov h , M. Hayek e ' 4 , I. A. Kachaev d , W. Kern 6 , 
E. King , 0. L. Kodolova f , V. L. Korotkikh f , M. A. Kostin f , 
J. Kuhn 1 , R. Lindenbusch a , V. V. Lipaev d , J. M. LoSecco h , 
J. J. Manak h ' 5 , J. Napolitano 1 , M. Nozar 1 , C. 01chanski c , 
A. I. Ostrovidov f , T. K. Pedlar §, A. V. Popov d , D. R. Rust a , 
D. I. Ryabchikov d , L. I. Sarycheva f , A. H. Sanjari h , E. Scott a , 
N. Shenhav e ' 4 , W. D. Shephard h , M. R. Shepherd a , 
N. B. Sinev f , K. K. Seth§, J. A. Smith 1 , P. T. Smith a , 
D. L. Stienike h , T. Sulanke a , A. Szczepaniak b , S. A. Taegar h ' 6 , 
D. R. Thompson 11 , I. N. Vardanyan f , D. P. Weygand c ' 5 , 
D. B. White 1 , H. J. Willutzki , J. Wise§, M. Witkowski 1 , 
A. A. Yershov f , and D. Zhao g 

^Department of Physics, Indiana University, Bloomington IN 47405, USA 

b Nuclear Theory Center, Indiana University, Bloomington IN 47405, USA 

c Department of Physics, Brookhaven National Laboratory, Upton, L.I., NY 11973 

^Institute for High Energy Physics, Protovino, Russian Federation 

^Department of Physics, University of Massachusetts Dartmouth, North 
Dartmouth MA 02747, USA 

f Institute for Nuclear Physics, Moscow State University, Moscow, Russian 

Federation 

^Department of Physics, Northwestern University, Evanston IL 60208, USA 
^Department of Physics, University of Notre Dame, Notre Dame IN 46556, USA 
department of Physics, Rensselaer Polytechnic Institute, Troy NY 12180, USA 



Abstract 

The parameter a = —0.005 ± 0.007 (stat) ± 0.004 (syst) describing the density 
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of events within the r] — ► 3ir Dalitz plot has been measured using data from 
Brookhaven experiment E852. The result is compared to previous values and theo- 
retical predictions. 



1 Introduction 



Chiral perturbation theory offers a consistent description of low energy QCD 
particularly applicable to meson decays. An outstanding difficulty is the large 
discrepancy between the predicted and observed rate of the r/ ^ 3n de- 
cay (|T]). We report here a measurement of the transition matrix element 
| M(j] —+ 3ir°) | 2 . The density of events on the Dalitz plot has been predicted to 
be uniform by Gasser and Leutwyler (HD and by Di Vecchia (0). More recently 
a calculation of Kambor, Wiesendanger and Wyler (|3]) found a density which 
decreases as the distance from the center of the plot increases. 

The Dalitz plot density is specified by a single parameter, conventionally de- 
noted as a. Kambor et al. show that both a and the rate T(t] — > 3tt) are 
sensitive to a parameter of the theory ( c in their notation ). The dependence 
of a and T on c is such that a more negative value of a reduces the value of T 
thereby increasing the disagreement with experiment. Two values of a are de- 
termined by the authors of (|3]), -0.014 and -0.007, the latter containing higher 
order corrections and leading to a prediction for T(t] — > 3ir) more consistent 
with experiment. 

To second order in the center of mass pion energy (f|) - (|7|) the Dalitz plot 
density for the decay r\ — > 3n° is parameterized as 

\M\ 2 ocl + 2az (1) 
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with 




where Ei is the energy of the z'th pion in the r] rest frame, p is the distance 
from the center of the Dalitz plot and p m ax the largest kinematically allowable 
value of p. 

A uniformly populated Dalitz plot, (as predicted by (]l|) and (||)) gives a = 0. a 
has been measured previously by Baglin et al. (|5|) , Aide et al. (§) , and Abele et 
al. ([?]). Recent results are shown in figure 1 along with the result of this work. 
Reference ([!]) finds a value significantly different from zero and in probable 
disagreement (2a) with theoretical expectations. The results reported in (|5]) 
and (|^) are consistent with both zero and the non-zero expectation of @. 
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Fig. 1. Previous measurements of a and the result of this work with a comparison 
to recent theoretical predictions. Not shown is the result of reference (||). 

In this work we report on a determination of a from the decay rj — > 3tc° where 
the rj was produced by 18.3 GeV/c negative pions incident on liquid hydro- 
gen. The data were collected in 1995 by the Brookhaven National Laboratory 
E852 collaboration. The apparatus has been described previously and 
consisted of a large, segmented lead glass calorimeter (LGD), charged particle 
tracking, veto and trigger systems capable of determining the charged particle 
multiplicity and the total energy deposited in the lead glass. 
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The data used for this analysis were collected with the "all-neutral" trigger 
which required no charged particles, an energy deposition greater than 12 GeV 
in the lead glass and no photons in the downstream photon veto. 

T) — > 37T° decays were selected by requiring exactly 6 reconstructed clusters in 
the LGD and a visible total energy greater than 16.5 GeV. All combinations of 
assignments to the hypothesis 3tt° were tested and if the best x 2 corresponded 
to a confidence level greater than 1% the event was considered a candidate. 
A preliminary fit using only the 7r° mass constraint was performed to allow 
evaluation of the 3ir° effective mass. Figure 2 shows the resulting distribution. 
Events with a mass less than 0.65 GeV/c 2 (corresponding to 3cr from the t] 
mass) were subjected to a full kinematic fit with the hypothesis n~p — > rjn; 
rj — > 37T°. A confidence level (CL) cut selected decays for the final analysis. 

It also was required that no photon have an energy less than 0.25 GeV and 
that the minimum separation between photons be larger than 9 cm. The effect 
and purpose of these requirements is discussed below. The final resulting data 
set contained 87,500 events. 

To determine a, eq. 1 is fitted to the observed distribution of z corrected 
for acceptance and phase space dependence. The acceptance correction was 
based on Monte-Carlo simulation, taking into account the apertures of the 
apparatus and a large sample of GEANT generated electromagnetic showers. 
The requirements on the minimum photon energy and the minimum allow- 
able photon separation mentioned above removed any uncertainties in the 
Monte-Carlo simulation associated with the transverse development of elec- 
tromagnetic showers and the behavior of low energy photons. The separation 
requirement also selected events where the nearest two electromagnetic show- 
ers were well resolved. The phase space dependence of z was removed by divid- 
ing the observed distribution by the distribution due to a uniformly populated 
Dalitz plot. 

The distribution resulting from removing the phase space dependence and 
correcting for the acceptance is proportional to M 2 and is shown in figure 3 
along with the result of the fit to eq. 1. The distribution has been numerically 
scaled to give a proportionality constant of one. Our result is a = — 0.0047 ± 
0.0074, where the error is statistical only. 

Systematic effects have been considered and the important contributors iden- 
tified. Since a is a static property of the 77 a measurement of a cannot depend 
on any data selection requirement, for example, momentum-transfer, separa- 
tion of photons in the detector, photon energy detection thresholds, confidence 
level of fits, etc. It was observed that our measured a remained constant as a 
function of the Confidence Level (CL) cut, down to about CL = 30%, below 
which a began to fall slightly. This is consistent with an increasing contribu- 
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Fig. 2. The observed 3tt° effective mass distribution. The r\ signal dominates the data 
with a secondary peak at ss 0.9 GeV/c 2 due to the decay K*° — > K°7r;K° — > 27r° 
visible. The inset shows the 77 region. There are 161,000 events under the peak 
with an estimated signal to noise ratio of 100 at the mass of the rj. The estimated 
background is indicated by the dashed line. 

tion from reconstructed events with one or more wrong- assignments of a 7 to 
a 7T°. It was required the an event have a confidence level greater than 30% for 
this analysis. Studies of the variation of a with the value chosen for the CL 
cut yielded an estimate of ±0.003 for the systematic contribution associated 
with this cut. 

The analysis further required that no two photons have reconstructed impact 
positions closer than 9 cm, a distance denoted by Ar. Variation of the value of 
Ar used in the analysis between 8.5 and 11 cm yielded our estimate of ±0.002 
for the contribution to the systematic error due to this effect. 

The effect of requiring the lowest energy photon to have an energy larger 
than 0.25 GeV was investigated by removing this requirement (increasing the 
data set by 1,140 events) and by replacing it with a requirement of 0.5 GeV 
(decreasing the data set by 15,230 events). A contribution of ±0.001 to the 
systematic error was found. 
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Fig. 3. Solid circles: z-distribution of 7] — > 37r° events on the Dalitz plot. The solid 
line is the result of a fit corresponding to eqn. 1. Open squares: the experimental 
acceptance as a function of z. 



For 20% of the events in the data sample it was possible to choose more 
than one assignment of photons to pions. The analysis presented here chose 
the assignment with the best confidence level. The effect of this choice was 
studied by choosing the second best assignment (when possible) and by using 
all assignments with an appropriate weight to calculate a. No statistically 
significant differences were observed. 



Finally, the data set was divided into two statistically independent halves by 
experimental run numbers and, separately, by range of momentum transfers. 
In both cases, the values determined for a were not statistically significantly 
different leading us to conclude these selections do not contribute to our sys- 
tematic error. All effects considered combine to give our final estimate of 
±0.004 for the systematic error. 
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2 Conclusion 



We have measured the a parameter of the transition matrix element M(i] — > 
3tt ) to be a = -0.005 ± 0.007 (stat) ± 0.004 (syst). Our result differs from 
one of the previous measurements (|?p by slightly more than 2a but is consis- 
tent with several other (|6]), ([5]) measurements. When compared to theoretical 
expectations, the value of zero favored by Di Vecchia, Gasser and Leutwyler 
([l|) , @, cannot be ruled out. The result reported here is consistent with the 
two values (-0.007 and -0.014) reported by Kambor, Wiesendanger and Wyler 

(D- 
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